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ABSTRACT 



We underscore the physics advantage of an experiment where neutrinos produced in a beta- 
X \ beam facility at CERN are observed in a large magnetized iron calorimeter (ICAL) at the India- 
c3 ' based Neutrino Observatory (INO). The CERN-INO distance is close to the so-called "magic" 
baseline which helps evade some of the parameter degeneracies and allows for a better measurement 
of the neutrino mass hierarchy and #13. We expound the possibility of using radioactive 8 B and 
8 Li as the source isotopes for the v e and v e beta-beam, respectively, and show that very good 
sensitivity to both the mass hierarchy and #13 is possible with a boost 7 in the 250-500 ballpark. 
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1 Introduction 



Spectacular results from a series of experiments involving solar [1], atmospheric [2], reactor [3], and 
accelerator [4, 5] neutrinos have firmly established neutrino oscillations and heralded the precision 
era in neutrino physics. Next generation experiments have been planned/proposed world-wide to 
further pin down the values of the oscillation parameters, Am^ 1 and 9± 2 for solar neutrinos and 
Am^ and #23 for atmospheric neutrinos. However, a vital task to be undertaken in the immediate 
future is the determination of the hitherto unknown mixing angle 9\s 2 . Discovery of a non-zero 
value for # 13 would open up the possibility of observing CP-violation in the lepton sector. Another 
outstanding problem concerns the ordering of the neutrino mass states, aka, the neutrino mass 
hierarchy. If Nature entertains a large enough mixing angle 13 , it could be possible to ascertain 
S(?n(Am| 1 ) and hence the neutrino mass hierarchy 3 . 

The neutrino mass hierarchy can be probed in terrestrial experiments through matter effects 
[8, 9, 10]. Large matter effects are known to exist in the oscillation probability v e — > [P e p)- 
Since it also provides the best way of determining the mixing angle #13 and the CP phase, Sep, 
it has been popularly hailed as the "golden channel" [11]. The P e ^ channel can be studied in 
experiments which use an initial u e (or u e ) beam and a detector which can efficiently see muons 4 . 
An absolutely pure and intense v e (or 9 e ) flux can be produced using "Beta Beams" (/3-beam) 
[12]. There is a possibility of such a facility coming up at the CERN accelerator complex and an 
international R&D effort is underway to check the feasibility of this proposal [13]. The Tevatron 
at FNAL has also been discussed as a potential accelerator for producing /5-beams. Most studies 
involving beta-beams have concentrated on using megaton water Cerenkov detectors for observing 
oscillations through muons produced via (or The reason is that these analyses consider 
low to medium values for the acceleration of the /3 unstable ions which create a relatively low 
energy (anti) neutrino beam. The detector therefore should be capable of observing the resultant 
low energy muons efficiently and unambiguously. Water Cerenkov detectors are particularly well 
suited for this. However, in contrast to most of the previous work involving /5-beams, which 
have mainly shown the physics potential of these kinds of experiments in determining $13 and 
5c p, in this paper we are especially interested in deciphering the neutrino mass hierarchy and #13 
simultaneously. To that end, we suggest a /3-beam with (anti)neutrino energies where we expect 
the largest earth matter effects, that is, energies in the multi-GeV regime. These high energy 
(or Vfj) can be very efficiently observed in large magnetized iron detectors, such as the envisaged 
ICAL detector at the India-based Neutrino Observatory [14]. 

A serious complication with all long baseline experiments involving the golden P eti channel 
arises from degeneracies which manifest in three forms: 

1 We define Am? = mf — m3. 

2 Currently, only an upper bound (sin 26 13 < 0.17 at 3er [6, 7]) exists. 

3 The neutrino mass hierarchy is termed "normal" ("inverted") if Am^ is positive (negative). Note that the 
terms normal and inverted "hierarchy" in this article refers to the neutrino mass ordering only, i.e., whether v 3 is 
heavier or lighter than 1/2, and our discussions are valid even for a quasi-degenerate neutrino mass spectrum. 

4 This is complementary to the standard accelerator beams where the initial (anti)neutrino flux consists of 
(or Vfj) and where the relevant channel is the probability P Me . Produced from decay of accelerated pions, these 
conventional (anti) neutrino beams suffer from an additional hurdle of an intrinsic v e (v e ) contamination, which 
poses a serious problem of backgrounds. A /3-beam is comprised of pure v e (or D e ). 
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• the (0i3, 5cp) intrinsic degeneracy [15], 

• the (s(/n(Am| 1 ), 5cp) degeneracy [16], 

• the (6*23, 7r/2 - 6*23 ) degeneracy [17]. 

This leads to an eight-fold degeneracy [18], with several spurious or "clone" solutions in addition 
to the true one and severely deteriorates the sensitivity of any experiment. It has been shown 
[18, 19] that the problem of clone solutions due the first two types of degeneracies can be evaded 
by choosing the baseline of the experiment equal to the characteristic refraction length due to the 
matter inside earth [8, 18, 19, 20]. This special value goes by the colloquial name "magic baseline" 
[19]. As we will discuss in detail later, at this baseline the sensitivity to the mass hierarchy and, 
more importantly, #i 3 , goes up significantly [19, 21], while the sensitivity to Sep is absent. 

Interestingly, the CERN-INO distance of 7152 km happens to be tantalizingly close to the 
magic baseline. This large baseline also enhances the matter effect and requires traversal through 
denser regions of the earth. Thus, for neutrinos (antineutrinos) with energies in the range 3-8 
GeV sizable matter effects are induced if the mass hierarchy is normal (inverted). A unique aspect 
of this set-up is the possibility of observing near-resonant matter effects in the v e — > channel. 
In fact, to our knowledge, what we propose here is the only experimental situation where near- 
resonant matter effects can be effectively used in a long baseline experiment to study the neutrino 
mass matrix. We show in this paper that the presence of this near-maximal earth matter effect 
not only maximizes the sensitivity to the neutrino mass hierarchy, it also gives the experiment 
an edge in the determination of the mixing angle 6*13. The increase in the probability P eti due 
to near-resonant matter effects, compensates for the fall in the /5-beam flux due to the very long 
baseline, so that one can achieve sensitivity to #13 and mass hierarchy which is comparable, even 
better, than most other proposed experimental set-ups. Therefore a /5-beam experiment with its 
source at CERN and the detector at INO could emerge as a powerful tool for a simultaneous 
determination of the neutrino mass hierarchy and #13. 

In [22] such an experimental set-up was considered for the first time and the physics potential 
explored. It was demonstrated that both the neutrino mass hierarchy as well as 6*13 may be probed 
through such a set-up. The ions considered for the /9-beam in that work were the most commonly 
used 18 A^e (for v e ) and 6 He (for v e ). One requires very high values of the Lorentz boost for these 
ions (7 ~ 10 3 ) because the energy E of the beam has to be in the few GeV range to enable 
detection in the ICAL detector, which is expected to have a threshold of about 1 GeV 5 . Such 
high values of 7, although possible in principle, might turn out to be very difficult to realize. 
Very recently two other ions, S B (for u e ) and 8 Li (for z/ e ), have been projected as viable options 
for the /3-beam source [23]. The main advantage that these ions offer is their substantially higher 
end-point energy, E . This allows one to access E ~ few GeV very easily with medium values of 7, 
that could be possible to achieve with either the existing CERN technology, or with the projected 
upgrades. The physics potential of these alternative ions to study standard neutrino oscillation 
parameters in the context of a CERN to Gran Sasso /3-beam experiment has been expounded in 
[24], while in [25] they have been studied to test R-parity violating SUSY. 

5 Most studies on the physics reach of /3-beam experiments have energies in the few hundred MeV range. 
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Figure 1: The unoscillated v e spectrum at INO-ICAL from the 8 B (left-hand panel) and 18 iVe 
(right-hand panel) /9-beam source at CERN, for three different choices of the boost factor 7 = 250 
(black dot-dashed curve), 350 (blue dashed curve), and 500 (red solid curve). 

In this paper we will work with 8 B and 8 Li as the candidate ions for the /3-beam source at 
CERN and ICAL at INO as the far detector. We perform a \ 2 analysis of the projected data in 
a future CERN-INO /3-beam experiment and present our results for three plausible benchmark 
values of the Lorentz boost factor 7. We begin in section 2 with a brief description of the /3- 
beam produced at the source. In section 3 we discuss the oscillation probability and highlight the 
importance of the magic baseline and the near-resonant matter effects in the P e/1 channel. Section 
4 contains our results for the event rate expected in the ICAL detector. In section 5 we present 
the details of the statistical analysis used. The sensitivity of the results on the mixing angle #13 
are presented in section 6, while section 7 focuses on the sensitivity to the neutrino mass hierarchy. 
In section 8 we consider exploiting the potential of muon charge identification at ICAL at INO 
through a /5-beam running simultaneously in both the u e and u e mode. We end in section 9 with 
our conclusions. 

2 The /3-Beam Fluxes 

The concept of a /?-beam was proposed by Zucchelli [12, 26]. The idea is to produce a pure, intense, 
collimated beam of u e or v e through the beta decay of highly accelerated and completely ionized 
radioactive ions circulating in a storage ring [27]. The resultant beam therefore comprises only 
of v e (or i? e ) and has definite characteristics, determined by the well known beta decay spectrum 
and the Lorentz boost factor. The former feature results in very low backgrounds while the latter 
ensures that there are essentially negligible flux systematics. This proposal is being studied in 
depth and fully exploits the existing CERN accelerator complex. The main future challenge lies 
in manufacturing an intense proton driver and a hippodrome-shaped decay ring which are at the 
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Figure 2: Same as in Fig. 1, but for the Vf, spectrum from 8 Li (left-hand panel) and 6 He (right- 
hand panel) source. 

core of this program. 



2.1 The Beta Decay Spectrum 

In the rest frame of the radioactive ion the beta decay spectrum is given by 

d 2 T* _ 1 
dQ*dE* ~ Anmlfr 



(E Q -E*)E*J(E -E*f-ml, 



where m e is the electron mass, Eq the electron total end-point energy, E* and r are the neutrino 
energy and lifetime of the decaying ion respectively in the latter's rest frame and / is a function 
of m e /E Q . In the lab frame, the flux of the unoscillated /3-beam at the detector is given by 



<j>(E, 



Na 



1 



47rL 2 ml f 7(1 — (3 cos ( 



(E -E*)E* 2 J(E -E*y 



rat 



(2) 



where L is the distance between the source and detector, Np are the number of useful decays in the 
storage ring per unit time, 9 is the angle between the neutrino flight direction and the direction 
in which the ions are boosted 6 and 7 is the Lorentz boost such that E* = 7-E(l — /3cos6), E 
being the neutrino energy in the lab frame. The maximum energy of neutrinos produced by a 
/3-beam with a Lorentz factor 7 is given by 



E u 



(E - m e ) 
7(1 -(3cos0)' 



(3) 
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Ion 


r(s) 


E (MeV) 


/ 


Decay fraction 


Beam 


18 Ne 


2.41 


3.92 


820.37 


92.1% 


v e 


tRe 


1.17 


4.02 


934.53 


100% 


Ve 


1* 


1.11 


14.43 


600684.26 


100% 


Ve 


|Li 


1.20 


13.47 


425355.16 


100% 





Table 1: Beta decay parameters: lifetime r, electron total end-point energy E Q , /-value and decay 
fraction for various ions [29]. 

2.2 Candidate Ions for the /3-beam 

From Eq. (2) it is seen that the total flux and energy of the /5-beam at the far detector depends 
mainly on the end point energy E of the beta decay ion and the Lorentz boost factor 7. The flux 
increases as both Eq and 7 increase. Larger 7 results in better collimation of the impinging flux, 
thereby increasing the statistics. The spectrum also shifts to larger energies as E and 7 increase. 
The neutrino cross section in the detector increases with energy, so for the same total flux a harder 
spectrum further enhances the statistics. Since the flux falls as 1/L 2 as the source-detector distance 
increases, high values of either E or 7 or both are needed to have sizable number of events at the 
far detector. E is an intrinsic property of the decaying ion while 7 is restricted by the nature of 
the accelerators and the machine design. The other desired properties which the ion should have 
include high production yield, large decay fraction, and a life time that is long enough to allow 
the ions to be adequately accelerated. It is also easier to store larger number of lower- Z isotopes 
in the storage ring [28]. In Table 1 we show the characteristic features of the four different ions 
which have been discussed in the literature as possible candidates for the /9-beam. While 18 Ne and 
8 B are (3 + emitters (producing a v e beam), 6 He and 8 Li are f3~ emitters (producing a v e beam). 
They all have comparable life times, conducive to the requirements necessary for the /3-beam, very 
high (or maximal) decay fraction and very low A/Z ratio. Note that the end-point energies, Eq, 
for 8 B and 8 Li are much larger than those for 18 Ne and 6 He. 

In this paper we are interested in the physics potential of an experimental set-up with a (5- 
beam at CERN and a large magnetized iron detector in India. For such a large baseline, one needs 
a very high value of 7, both to cross the detector energy threshold as well as to get reasonable 
statistics in the detector. In particular, with 6 He (or 18 iVe) as the source ion, one needs 7 ^ 1000. 
Such high values of 7 can only be achieved by using the LHC itself. As noted above, E for 
the alternative options for the radioactive ion source ( 8 B and 8 Li) is higher by a factor of about 
four. This means that it should be possible to produce high intensity, high energy beams with 8 B 
and 8 Li for a much lower boost factor 7 . We show in Fig. 1 the 8 B (left-hand panel) and 18 Ne 
(right-hand panel) /3-beam flux expected at INO as a function of the v e energy, for three different 
benchmark values of 7. Fig. 2 shows the corresponding spectra for the 8 Li (left-hand panel) and 
6 He (right-hand panel) v e /3-beam. Note that even though apparently it might seem from the 

6 We work with the on-axis flux for which 9 = 0. 

7 The loss in collimation is not significant as we show later. 
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= 2.5 x 10" 


3 eV 2 


sin 2 29 23 


= 1.0 




Amli = 


8.0 x 10~ 5 


eV 2 


sin 2 12 = 


= 0.31 




5cp = 







Table 2: Chosen benchmark values of oscillation parameters, except sin 2 26 ) i 3 . 



figures that the v e (p e ) flux is larger for 18 Ne ( 6 He), in reality, for a given 7, the total flux is given 
by the area under the respective curves. One can easily check that for a fixed 7 this is same for 
both the ions as expected, since we have assumed equal number of decays for both 18 Ne and 8 B 
for v e ( 6 He and 8 Li for z/ e ) 8 . Also note that even though the total v e (P e ) flux remains the same for 
both the ions for a given 7, the number of events produced in the detector is much enhanced for 
the 8 B ( s Li) beam since the energy of the beam is larger and the charged current cross sections 
increase with the neutrino energy. Studies have shown that it is possible to accelerate 6 He to 
7 ^ 250 with the existing facilities at CERN, while 7 = 250 — 600 should be accessible with 
the "Super-SPS", an upgraded version of the SPS with super- conducting magnets [30, 31]. The 
Tevatron at FNAL could in principle also be used to produce a /5-beam with 7 < 600. 

In the low 7 design of beta-beams 2.9 x 10 18 6 He and 1.1 x 10 18 18 Ne useful decays per year 
should be possible to achieve. Earlier, only these ions were considered because it was believed that 
8 B could not be produced with the standard ISOLDE techniques. Since most exercises focused on 
observation of CP-violation, it was necessary to have both v e and v e beams with similar spectra, 
so 8 Li (though considered in [28]) was also generally ignored. Interest in both these ions have been 
rekindled recently [23, 24], as it appears that having intense S B and 8 Li fluxes should be feasible 
using the ionization cooling technique [23]. In what follows, we will vary the value of 7 to test 
the physics potential of the CERN /3-beam INO-ICAL set-up and, following the current practice, 
assume that it is possible to get 2.9 x 10 18 useful decays per year for 8 Li and 1.1 x 10 18 for 8 B for 
all values of 7. 



3 Oscillation Phenomenology — The Magic Baseline 



In this section we will very briefly review some issues related to the neutrino oscillation probability. 
We consider the standard form of the neutrino mixing matrix 

/ C12C13 S12C13 s 13 e~ l5cp \ 



U 



C12C13 

-S12C23 - Ci 2 S 2 3Sl3e 



iScr 



V S12S23 - Ci2C 23 Sl3e 



iSci 



S12C13 

C12C23 - Sl2S23Sl3e 
-C12S23 - Sl2C 2 3Sl 3 e 



iScp 
iScr 



S23C13 
C23C13 J 



(4) 



3 Of course, the on-axis flux increases with 7 because of better collimation of the beam. 
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Figure 3: The transition probability P efJL as a function of E for four values of the baseline L. 
The band reflects the effect of the unknown 5c p- The dark (red) shaded band is for the normal 
hierarchy (NH) while the light (cyan) shaded band is for the inverted hierarchy (IH). We have 
taken sin 2 2^13 = 0.1 and for all other oscillation parameters we assume the benchmark values 
given in Table 2. 



where c^- = cos Sij = sin0y and we ignore the Majorana phases. From the analyses of the 
current solar, atmospheric, reactor, and accelerator data set one has the following 3a constraints 
on the oscillation parameters [7]: 



7.2 x l(T 5 eV 2 < Am 2 x < 9.2 x l(r 5 eV 2 and 0.25 < sin 2 # 12 < 0.39 



2.0 x 10~ 3 eV 2 < | Am 2 x | < 3.2 x 10" 3 eV 2 and sin 2 29 23 > 0.9 



sin 2 2^13 < 0.17 



(5) 
(6) 
(7) 



Throughout this paper we show all our results assuming certain true values for the oscillation 
parameters lAm^l, sin 2 #23, Am^, sin 2 9 12 and 5cp given in Table 2. We will let sin 2 2^13 take 
all possible true values (which we will call sin 2 20 13 (true)) and present our result as a function of 
sin 2 26 , 13 (true). 
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Figure 4: The dark (red) shaded band is the same as in Fig. 3. The light (green) shaded band 
shows the corresponding P efl for sin 2 26?i3 = 0.05. Values of all the other oscillation parameters are 
same as in Fig. 3 and the hierarchy is assumed to be normal. 



Since both sin 2 2# 13 and a = Am^/Amfi are very small parameters, as can be seen from the 
current constraints, we can expand the so-called "golden channel" probability in the constant 
matter density approximation as 9 [11, 32] 

,„ o „ sin 2 [(l-i)Al 

P eu ~ sin 2 # 23 sin 2 28 13 l ± 

* {I -Ay 

sin(iA) sin[(l - A) A] 



+ a sin 26*13 sin 29\2 sin 26*23 sin 5cp sin(A) 



+ a sin 26*13 sin 26*12 sin 26> 2 3 cos 5cp cos(A) 

9 9 „ 9 n sin 2 (AA) 

+ a 2 cos 2 #23 sin 2 26*12 ^ 

^42 



A (1-A) 
sin(iA) sin[(l - A) A] 
A (1-A) 



9 This particular low order expansion of the transition probability is valid only in the range of L and E where the 
resonance condition A = 1 is never reached. For the E and L range that we consider in this paper, this condition 
is satisfied and Eq. (8) fails, as we discuss later in this section. Nonetheless, we present this expression to illustrate 
the effect of the magic baseline. 
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where 



±A 

l 31 



A = 2V2G F N e E being the matter potential, where N e is the electron number density inside the 
earth and we have expanded the probability in a and sin2# 13 , keeping only up to second order 
terms in both these parameters. The + ( — ) sign in Eq. (10) refers to neutrinos (antineutrinos). 
The first term of Eq. (8) can be used to extract information about the value of 6>i 3 . This is also the 
term which has the largest earth effect and this effect of matter can be used to determine the sign 
of Amjj. The second (CP violating) and third (CP conserving) terms depend on the CP phase 
5c p and can be used to find it. The last term is independent of both #13 and 5c p and depends 
mainly on the solar parameters Am^ and 12 . In very long baseline experiments, this term has 
sizable matter effects and if the true value of #13 turns out to be zero (or nearly zero), this would 
be the only surviving term in P eM which could still be used to study matter enhanced oscillations 
[33]. 

Despite its advantage in determining the most interesting oscillation parameters, this channel 
is, however, rife with the problem of "degeneracies" - the 9\^-5cp intrinsic degeneracy [15], the 
sgn^Am^) degeneracy [16] and the octant of 6*23 degeneracy [17] - leading to an overall eight- 
fold degeneracy in the parameter values [18]. Various schemes to tackle some or all of these 
problems have been studied in great detail in the literature. These can be broadly categorized 
into ones where data from experiments at different energies and different baselines are combined 
together [15, 34, 35, 36, 37, 38], where data from accelerator experiments observing different 
oscillation channels are considered together [39, 40], and where data from atmospheric [41, 42] or 
reactor experiments [43] are combined with the golden channel measurements in an accelerator 
experiment. 

A particularly interesting scenario arises when the condition 

sin(iA) = (11) 

is satisfied. In such an event, the last three terms in Eq. (8) drop out and the P eil channel enables 
a clean determination of #13 and s^n(Am| 1 ). 

Since A A = ±AL/AE by definition, the condition (11) reduces to pL = y/2ir/G F Y e , where Y e 
is the electron fraction inside the earth. This gives 

1 ~ 32725 , (12) 



[g/cc] [km 

which for the PREM [44] density profile of the earth is satisfied for the "magic baseline" [19] 

L magic ~ 7690 km. (13) 
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We refer the reader to [20] for a very recent and enlightening discussion on the physical meaning 
of the magic baseline. 

The CERN-INO distance corresponds to L — 7152 km, which is tantalizingly close to the 
magic baseline. We therefore expect that the /3-beam experiment we consider here should give an 
essentially degeneracy-free measurement of both $13 and sg , n(Am| 1 ). 

The large CERN-INO baseline, of course, results in very significant earth matter effects in the 
P eM channel. In fact, for the baseline of 7152 km, the average earth matter density calculated 
using the PREM earth matter density profile is p av = 4.13 gm/cc, for which the resonance energy 

E res = |Am ' / i |cOS2gl3 (14) 
2y/2G F N e 1 ' 

= 6.1 GeV , (15) 

for [Amf-J = 2.5 x 1CT 3 eV 2 and sin 2 26 , i 3 = 0.1. One can check (cf. Figs 1 and 2) that this is 
roughly in the ballpark where we expect the maximum /9-beam flux for 8 B and 8 Li ions with the 
benchmark 7 values used in this paper. 

Note that the low order expansion of the probability P eil given by Eq. (8) is valid only for values 
of E and earth matter density p (and hence L) where flavor oscillations are far from resonance, 
i.e., A <C 1. In the limit A ~ 1, one can check that even though the analytic expression for P e/X 
given by Eq. (8) remains finite, the resultant probability obtained is incorrect [45]. We reiterate 
that Eq. (8) was presented only in order to elucidate the importance of the magic baseline. For all 
the numerical results presented in this paper, we calculate the exact three generation oscillation 
probability using the realistic PREM [44] profile for the earth matter density. 

The exact neutrino transition probability using the PREM density profile is given in Fig. 3 
as a function of the neutrino energy, for four different baselines (four panels). We allow Sep to 
take on all possible values between 0-360 degrees and the resultant probability is shown as a band, 
with the thickness of the band reflecting the effect of S C p on P efl . The figure is drawn assuming 
the benchmark values of the oscillation parameters given in Table 2 and sin 2 2# 13 = 0.1. We show 
the probability for both the normal hierarchy (dark band) as well as the inverted hierarchy (light 
band). As discussed in detail above, for L = 7500 km, which is close to the magic baseline, the 
effect of the CP phase is seen to be almost negligible, while for all other cases the impact of Sep 
on P eM is seen to be appreciable. In fact, for L = 1000 km, the probability corresponding to the 
normal and inverted hierarchy become almost indistinguishable due to the uncertainty arising from 
the unknown value of Sep. As the baseline is increased, earth matter density increases, enhancing 
the impact of matter effects. The probability for normal hierarchy is hugely enhanced for the 
neutrinos, while for the inverted hierarchy matter effects do not bring any significant change. This 
difference in the predicted probability, evident in the panels corresponding to L = 4000, 7500 and 
10000 km of Fig. 3, can be used to determine the neutrino mass hierarchy. 

In Fig. 4 we display the dependence of the neutrino probability P e/I on the mixing angle #13 for 
four different baselines. The dark bands, as in Fig. 3, are for normal hierarchy and sin 2 2#i3 = 0.1 
with full variation of Sep, while the light bands are for normal hierarchy and sin 2 29 13 = 0.05. The 
impact of matter effect in increasing the 6*13 sensitivity of a given experimental set-up is evident 
from the figure. The #13 sensitivity for L = 1000 km can be seen to be much weaker than for the 
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other cases, since matter effects are smaller. However, the most striking feature seen in Fig. 4 is 
the effect of the magic baseline in enhancing the sensitivity of the experiment to 9 13 . The figure 
clearly shows that the difference in the predicted P e)l for the two values of sin 2 26 1S is largest for 
L = 7500 km (since effect of 5qp is the least) and thus an experiment at this baseline is most 
suitable for probing #13. Figs. 3 and 4 therefore reinforce our choice of the near- magic baseline 
as one of the best options for determining the neutrino mass hierarchy and #13, since both these 
parameters are directly related to large matter effects and the uncertainty of Sep could prove to 
be a hindrance in their measurement at non-magic baselines. 

4 Event Rates in INO-ICAL 
4.1 The ICAL Detector at INO 

The proposed large magnetized iron calorimeter at the India-based Neutrino Observatory [14] is 
planned to have a total mass of 50 kton at startup, which might be later upgraded to 100 kton. 
The INO facility is expected to come up at PUSHEP (lat. North 11.5°, long. East 76.6°), situated 
close to Bangalore in southern India. This constitutes a baseline of 7152 km from CERN. The 
ICAL detector will have a modular structure with a total lateral size of 48m x 16m, divided into 
three modules of 16m x 16m each. Each of these modules will have 140 horizontal layers of ~ 6 
cm thick iron plates, separated from each other by a gap of ~ 2.5 cm to hold the active detector 
material, giving a total height of 12 m for the full detector. The active detector elements will 
be resistive plate chambers (RPC), made from either glass or Bakelite and will be filled with a 
suitable gas mixture, which will be recycled with approximately one volume change per day. An 
external magnetic field of ~ 1.3 Tesla would be applied over the entire detector. The detector will 
be surrounded by an external layer of scintillator or proportional gas counters which will act both 
as veto to identify external muon backgrounds as well as to identify partially contained events. 

According to the detector simulation performed by the INO collaboration, the detector energy 
threshold for jj^ is expected to be around ~ 1 GeV and charge identification efficiency will be 
about 95%. In what follows, we will present our numerical results assuming an energy threshold 
of 1.5 GeV, detector charge identification efficiency as 95% and unless stated otherwise, detection 
efficiency as 60% for fi^ (cf. Table 3). The 60% detection efficiency we have used is an extremely 
conservative estimate and it should be possible to increase it to at least 80%. We shall therefore 
also present our main results for this more optimistic estimate for the detector efficiency. However, 
unless it is mentioned otherwise, all our results correspond to the detector characteristics given in 
Table 3. We have explicitly checked that our results remain unaffected if the energy threshold is 
raised to 2 GeV for the entire range of assumed Lorentz boost factor 7 = 250 — 500. For 7 > 350 
the threshold can be even 3 GeV, while for 7 > 500 one can work with an energy threshold of 4 
GeV, without changing the final results. This can be seen from Figs. 1 and 2; for 7 = 350(500), 
the majority of neutrinos arriving at INO-ICAL would have E > 3(4) GeV. 
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Total Mass 


50 kton 


Energy threshold 


1.5 GeV 


Detection Efficiency (e) 


60% 


Charge Identification Efficiency (fin) 


95% 



Table 3: Detector characteristics used in the simulations. 



4.2 Oscillation Signal 

The total number of events collected in INO-ICAL from a u e (or u e ) /3-beam exposure over 

a period of T years is given by 

£max 

N^ = Tn n f ID e / dE dE A (/>(E) a^(E) R(E, E A ) P e/X (E) (16) 

where n n are the number of target nucleons, fm is the charge identification efficiency, e is 
the detection efficiency, 4>(E) (given by Eq. (2)) is the /3-beam flux at INO-ICAL in units of 
m _2 year~ 1 GeV _1 , a V)l is the detection cross section for z/ M in units of m 2 , E is the true energy 
of the incoming neutrino and Ea is the measured energy 10 ' 11 , R(E,Ea) is the detector energy 
resolution function 12 , _E max is the maximum energy of the neutrinos for a given Lorentz factor 7 
(given by Eq. (3)) and P efl is the v e — > z/ M oscillation probability. The expression for the /z + signal 
in the detector from a v e /3-beam flux is given by replacing <p Ue by 0p e , a Vfi by <Tp M and P efJi by P S p,. 
For the neutrino-nucleon interaction we consider quasi-elastic scattering, single-pion production 
and deep inelastic scattering and use the cross sections given in the Globes package [46] which are 
taken from [47, 48]. 

We show in Fig. 5 the number of events expected in INO-ICAL from 5 years exposure of 
an 8 B (for u e , shown in the left-hand panel) or s Li (for z/ e , shown in the right-hand panel) f3- 
beam from CERN. The expected number of events is presented as a function of sin 2 2#i 3 for 
both normal (NH) and inverted (IH) hierarchy for three benchmark values of the boost factor 7. 
Large resonant matter effects in the neutrino channel for normal hierarchy drives the number of 
expected events to very large values, compared to what would be expected for inverted hierarchy. 
Similarly, in the antineutrino channel we have resonant matter effects for inverted hierarchy and 
the number of predicted events is many times larger than for normal hierarchy 13 . This difference 
in the number of events is seen to increase with sin 2 26> 13 up to a certain (large) value, beyond 

10 The "measured" energy of the neutrino can be reconstructed from the total measured energy of the muon and 
the accompanying hadron shower. 

n The integration over the reconstructed (measured) neutrino energy is done from the threshold energy E A h 
(taken as 1.5 GeV in this paper) to infinity. 

12 We assume a Gaussian resolution function with a = 0.15-E". 

3 In fact, matter effects are seen to suppress the oscillation probability for neutrinos (antineutrinos) compared 
to that in vacuum when the hierarchy is inverted (normal) . 
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Figure 5: The expected number of events in 5 years running time, as a function of sin 2 2# 13 . The 
value of 7 and the hierarchy chosen corresponding to each curve is shown in the figure legend. 
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Figure 6: The iso-event curves in the 5cp-sin 2 2#i3 plane for four baselines are shown in the figure. 
The true values of Sep and sin 2 2#i3 are assumed to be 0° and 0.043 respectively. Left- (right-) 
hand panel is for the v e (u e ) /3-beam. The assumed hierarchy is mentioned in the figure. 
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which it starts to decrease. This limiting value of sin 2 2#i 3 at which the number of events in the 
neutrino (antineutrino) channel peaks for the normal (inverted) hierarchy comes from an interplay 
of two effects which we now discuss. In the approximation where Am^ can be neglected, the 
probability P efl would be largest if the mixing angle in matter (sin 2 29f^) and the mass squared 
difference driven oscillatory term in matter (sin 2 [(Am| 1 ) M L/4£ 1 ]) are maximum simultaneously. 
In other words, largest matter effects come when the resonance energy E res is almost equal to the 
energy at which the oscillatory factor goes to one [49, 50, 51]. For a given L (and hence matter 
density, p), one can calculate the value of sin 2 2# 13 for which this condition can be satisfied and it 
is given by [50, 51] 

16.18 x 10 3 

tan2# 13 = - — 17 

p(gm/cc)L(km) 

For L = 7152 km, it turns out that this value is sin 2 2#i 3 ~ 0.23. Therefore, matter effects 
and hence the event rate in the detector keeps increasing until this value of sin 2 26>i 3 is reached. 
Beyond this limiting value of sin 2 26*i 3 , effect of earth matter in P eM falls and the event rate for 
the neutrinos, for normal hierarchy, and antineutrinos, for inverted hierarchy, decreases. However, 
note that this value of sin 2 29 13 is already disfavored from the CHOOZ data [6, 7]. 

From the left-hand panel of Fig. 5 we note that for 7 = 500 and sin 2 2# 13 = 0.05, the predicted 
number of neutrino events for normal hierarchy is 101, while that for inverted hierarchy is only 3. 
This implies that if the normal hierarchy was true, we could comprehensively rule out the wrong 
inverted hierarchy. For smaller boost factors 7, the difference between the number of events 
decreases; however, we could still determine the neutrino mass hierarchy if sin 2 2#i 3 (true) is not 
very small. We will perform a statistical analysis of projected data in the following sections and 
present our results on the hierarchy sensitivity in section 7. Instead of the neutrino, one could 
work with the V e /3-beam and achieve similar sensitivity to the neutrino mass hierarchy. Note 
that while the interaction cross section for the z/ e s are much smaller, the flux itself is larger owing 
to the (assumed) larger number of decays per year for 8 Li. Thus the statistics expected in both 
the neutrino as well as the antineutrino channel is comparable and, as we will see, the hierarchy 
sensitivities for the v e and v e /5-beam are hence similar. 

As discussed above, resonant matter effects in the neutrino (antineutrino) channel for the 
normal (inverted) hierarchy, results in substantial enhancement in the observed number of events. 
In particular, we note from Fig. 5 that the number of events in the near-resonant channels depend 
strongly on the value of sin 2 2# 13 , since the extent of matter effects is dictated directly by #i 3 . We 
can see from the figure that the dependence of the event rate on sin 2 2# 13 is much enhanced due to 
matter effects. Therefore, if the true hierarchy was normal (inverted), we could use the neutrino 
(antineutrino) /3-beam in the proposed experiment to measure/constrain the mixing angle 6>i 3 . 

In Fig. 6 we show the effect of Sep on the variation of the measured rate with the value of 
the mixing angle #13, for four different baselines 732 km, 4000 km, 7152 km and 10000 km. The 
left-hand panel shows the results for the v e /3-beam assuming a normal mass hierarchy, while the 
right-hand panel gives the same for the v e /3-beam flux with inverted hierarchy. For all cases in 
any panel, we have considered a 50 kton magnetized iron calorimeter as the far detector and the 
same flux created at the source. Each curve gives the sets of values of {sin 2 2#i 3 , Sep} which give 
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the same observed rate in the detector as the set {sin 2 26*i3 = 0.04, 5cp = 0°}. In other words, 
if the true value of sin 2 29 13 and Sep were 0.04 and 0° respectively, then every point on a given 
curve would also be a solution for that experiment. It is clear from this figure that combining 
results from experiments at different baselines helps solve/reduce the problem. However, the most 
important issue exemplified here is the fact that for the baseline 7152 km, which is the CERN-INO 
distance, the effect of the unknown value of 5cp on the measurement of #13 and the mass hierarchy, 
is very small. This happens because this distance corresponds to a near-magic baseline for which, 
as noted earlier, the 5c p dependent terms are almost vanishing. 

4.3 Backgrounds 

The possible backgrounds in a v e /3-beam experiment 14 using /i~ as an oscillation signal come from 
neutral current events such as 

v x + d{u) — > v x + d(u) (18) 
v x + d(u) — > v x + d(u) + qq (19) 

and u e charged current events 

v e + d — > e~ + u or c (Cabibbo suppressed) (20) 
v g + d — > e~ + u + qq (21) 

The quarks in the final state could produce mesons as a part of the hadronic junk. These mesons 
may then decay producing secondary muons, giving rise to a signal which might constitute a 
possible background. 

As discussed earlier, INO-ICAL will have 6 cm thick iron plates. Such a dense tracking detector 
will have excellent muon/pion and muon/electron separation capability in the energy range we are 
working with. The simulations carried out by the INO collaboration have shown that after the 
standard kinematical cuts are imposed, the electrons do not give any signal at all, while in 99% 
of the cases, the pions and kaons get absorbed very quickly in the iron for the energy range of 
interest to us and therefore do not hit enough RPCs to give a signal. At the energies of the /?-beams 
considered here, production cross section of D mesons (also Cabibbo suppressed) is small and they 
do not constitute a problem for the experiment. The associated strange or charm production is 
also highly suppressed at these energies. In addition, the fact that the detector will have a charge 
identification capability means that secondary /i + produced can be safely discarded, reducing the 
background even further. Therefore, in this analysis we do not consider any backgrounds coming 
from either the neutral current events or charged current events of f e 15 . 

Since the oscillation probability v e — > v T is about the same as that for v e — > z/ M , we expect 
almost as many v T arriving at the detector as v^. The r~ produced through charged current 

14 The discussion concerning the v e /3-beam is similar and hence is not repeated here. 

15 Mesons produced in neutral current processes are degraded in energy. Note that backgrounds from these 
mesons are very important in the case of the neutrino factory. However, since our relevant energy range is lower, 
the mesons produced in each event are much lower in energy and hence can be easily rejected by putting suitable 
cuts. In our analysis we have assumed cuts that are stringent enough to completely reject these backgrounds. 
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interaction may decay producing secondary ^ with a branching ratio of 17.36%. But, the r 
threshold (3.5 GeV) is high and the production cross section suppressed compared to fi. So we do 
not expect any significant background from this source either. We have estimated the number of 
secondary muons produced from the v T component of the beam. For sin 2 2# 13 = 0.01, we expect 
0.008, 0.061 and 0.2 muon events per year respectively, for 7 = 250, 350 and 500. In addition, 
these secondary muons will be severely degraded in energy and therefore can be eliminated through 
energy cuts. We therefore neglect the backgrounds from this source as well. 



5 Details of the Statistical Method 

In order to quantify the sensitivity of this /5-beam experimental set-up to the mixing angle #13 
and sgn(Am^ 1 ), we perform a statistical analysis of the "data" generated in INO-ICAL, assuming 
certain true values of the parameters. Since the number of events expected in the detector might 
be very small, depending on how small the true value of sin 2 29i 3 is, we define a x 2 assuming 
Poissonian distribution for the error as, 



X 2 (M) = min 

Zk 



2 N 



th 



N ex In 



fjth 



)+E3 

/ k J 



(22) 



In Eq. (22), {00} is a set of oscillation parameters, N ex is the observed number of events, while 
the systematic errors in the data and the theory are accounted for through the set of "pulls" 
where k runs over the different systematic uncertainties involved. The pulls are defined in such a 
way that the number of expected events, N th , corresponds to £ k = 0. N th is the number of events 
when the systematic errors are included in a manner such that the effect of the fc-th uncertainty 
at the ±1(7 level is picked up when = ±1: 



N th ({u J },{^}) = N^{u}) 



thi 



K 



k=i 



(23) 



In Eq. (23) N th ({ou}, has been expanded in powers of keeping only linear terms. The 
quantities ir k give the fractional rate of change of N th ({uj}) due to the kth systematic error. In our 
analysis, we include systematic uncertainties in the normalization of the /3-beam flux at the source, 
the error in the cross section and the detector systematic uncertainty. For the flux normalization 
error we adopt a total uncertainty of 2%. The neutrino-nucleon cross sections, especially at large 
energies and for massive target nuclei, are known to have large uncertainties. We include a 10% 
error coming from the uncertainty in the interaction cross section. A total detector systematic 
uncertainty of 2% is also included. We do not include any systematic error related to the shape 
of the energy spectrum since we work only with the energy integrated total rates, for which such 
types of errors would be negligible. 

Unless otherwise stated, the data are always generated at the benchmark values of the oscil- 



lation parameters jAm^j 



sin 2# 23 , sin 9 12 and 5c p given in Table 2. The true value of 



9 13 as well as sgn^Am^) will be allowed to change, since these are the parameters under study in 
this experimental set-up of interest. Throughout the paper we will use the notation sin 2 2^13 (true) 
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as the true value of this parameter chosen by Nature, and sin 2 2#i 3 as the fitted value. The data, 
simulated with the (assumed) true values of the oscillation parameters, are fitted by calculating 
N th ({u}) for any set of values for the oscillation parameters. We next calculate the chi-square 
function x 2 (W}, {£fc}) f° r every possible N th ({cu}, {£*;})> obtained by varying the oscillation pa- 
rameters and We first minimize this x 2 ({u;}, {£&}) function with respect to the pulls This 
X 2 {{uj}) (given by Eq. (22)) is next minimized with respect to all the oscillation parameters 16 , 
which are permitted to assume all allowed values in calculating N th ({uj}). 

It is expected that the parameters |AmU and sin 2 26* 2 3 would be fairly well determined by the 
long baseline experiments T2K [53] and NOz/A [54]. In order to take into account this information, 
which should be available by the time the /3-beam facility comes up at CERN, we add to our x 2 
the "prior" function such that 

(24) 



2 / 1 Am§! | - | Am§! (true) | \ 2 / sin 2 26» 23 - sin 2 26» 23 (true) 
X ^or~y ' a(Aml) ) + [ a(sm 2 26 23 ) 

where Am 2 1 (true) and sin 2 2#23 (true) are taken as the benchmark values in Table 2. We will assume 
that the la error on these parameters would be reduced to a(Am^ 1 ) = 1.5% and a (sin 2 26*23) = 
1.0% [55]. We minimize the total function Xtotai with respect to the oscillation parameters and 
present the C.L. results assuming as a convention that no corresponds to 

A Xtotai = n 2 , (26) 

where Ax 2 ota i is the minimum Xtotai obtained in the analysis for a given value of the parameter 
sin 2 2^13 (true) and/or the mass hierarchy. 

Note that we have not taken into account any error coming from the uncertainty in the matter 
density. Fig. 7 shows our justification for this. This figure shows the expected event rate for the 
v e /9-beam as a function of sin 2 20 13 . The solid curve corresponds to the PREM profile, while the 
dashed and dot-dashed curves are for PREM densities scaled by -5% and +5% respectively. We 
can see that the effect of the 5% uncertainty in the matter density profile shows up in the event 
rate only for values of sin 2 2# 13 already disfavored by the CHOOZ experiment. In particular, we 
note that for sin 2 2# 13 ^0.1 the effect of the matter density uncertainty is negligible. 




6 Determining sin 2 2#i3 

We begin by presenting results on how well sin 2 2#i 3 can be determined by the set-up we are 
proposing. We divide our results into two subcategories. First we assume that sin 2 2#i 3 (true) = 

16 In order to save computer time, we keep Am^ and sin 2 9 12 fixed in our analysis. There is two-fold motivation for 
this approximation. Firstly, these parameters, especially Am 2 ,!, i s expected to be measured with a high precision in 
the proposed solar and reactor neutrino experiments [52]. Secondly since we are working at a near-magic baseline, 
we do not expect any significant impact of these parameters on our results. 
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sin 20i3 

Figure 7: Number of events versus sin 2 2#i3 for the v e /3-beam with 7 = 500, showing the impact of 
uncertainty in the earth matter density profile. The solid curve corresponds to the PREM matter 
density profile, the dashed (dot-dashed) curve is for the PREM density scaled by -5% (+5%). 



and expound the possibility of improving the upper bound on this parameter. This is what we 
would call the sensitivity of the experiment to sin 2 26 l i3. Next we will assume that sin 2 2# 13 (true) 
is large enough to give a positive signal in the experiment. In that case, we will probe quanti- 
tatively how precisely this parameter could then be determined. All results are presented after 
marginalizing over lAmJjJ, sin 2 2#23, and Sep, as described in section 5. 

6.1 Sensitivity to sin 2 2^13 

We define the sensitivity to #13 as the minimum value of sin 2 2# 13 which this experiment would be 
able to distinguish from sin 2 2^13 (true) = 0. The results of our analysis are presented in Fig. 8, 
where we show the upper bounds on sin 2 2#i3 that this experiment can impose at the 90% CL. 
(solid curve for 5 years running and dotted curve for 10 years running) and 3a CL. (long dashed 
curve for 5 years running and dashed-dotted curve for 10 years running). The left-hand panel is 
for a u e /3-beam assuming that normal hierarchy is the true hierarchy, while the right-hand panel 
is for an v e /3-beam with the assumption that the inverted hierarchy is true. We show our results 
as a function of the Lorentz boost 7. The 90% and 3a CL. upper bounds that we can impose on 
sin 2 2#i3 are displayed in Table 4 for the three benchmark values of 7 and with five and ten years 
of running of the experiment in the neutrino (columns 3 and 4) and in the antineutrino (columns 
5 and 6) mode. In the Table we show the sin 2 2# 13 sensitivity limits for the conservative detector 
efficiency of 60% for [V s1 as well as for 80% efficiency. This can be compared with the sensitivity we 
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Figure 8: The upper limit on sin 2 29\z that can be imposed by the CERN-INO /3-beam set-up 
with 5 years and 10 years of running. The left-hand panel shows the sin 2 2#i3 sensitivity expected 
for the neutrino mode assuming normal hierarchy to be true, while the right-hand panel shows 
the sensitivity for the antineutrino mode and an inverted hierarchy. 



expect from the other forthcoming/planned experiments. At the 90% C.L., the T2K and NOz/A 
experiments are expected to constrain sin 2 2#i3 < 2.3 x 10~ 2 and sin 2 2#i3 < 2.4 x 10 -2 respectively 
[55]. The Double-Chooz experiment is likely to push the limit down to sin 2 29 13 < 3.2 x 10~ 2 at 
90% C.L., which could be improved to sin 2 2# 13 < 9 x 10~ 3 by the Reactor-II set-up [55]. The 
"discovery potential" 17 of the CERN-MEMPHYS project for five years of running of both the f3- 
beam and the SPL super-beam set-ups considered alone is sin 2 26* 13 ^5x 10~ 3 at 3a [42]. With a 
combined data set in MEMPHYS for five years running in /3-beam and five years in the SPL super- 
beam, the discovery limit could be improved to sin 2 2#i3 ^ 3.5 x 10 -3 at 3cr [42]. The best limit, 
obtained after a thorough optimization of the /3-beam experimental options for the #13 sensitivity 
reach in [56] (see also [57]), is sin 2 29 13 ^ 1.5 x 10~ 3 at 3cr, for what the authors label as set-up 
l 18 and set-up 2 19 . Therefore, for the experimental set-up with the /3-beam source in CERN and 
INO-ICAL as the far detector the projected sin 2 29\ 3 sensitivity is comparable to, if not better, 
than most of the other planned schemes. In fact, the sensitivity reach of this experiment is even 
better than that expected for an entry level Neutrino Factory (NuFact-I) [58]. Only a NuFact-II 
with L = 7500 km [19] corresponding to the magic baseline, could give a sin 2 2# 13 sensitivity 

17 The authors define the "discovery potential" as the minimum value of sin 2 2#i3(true) for which sin 2 26>i3 = 
gives a A^ 2 = 9. 

18 This corresponds to successive 8 year runs of a /3-beam in the neutrino and antineutrino mode (total 16 years) 
with 7 = 200, L = 520 km and a water Cerenkov detector with 500 kton fiducial mass. 

19 This corresponds to successive 8 year runs of a /3-beam in the neutrino and antineutrino mode (total 16 years) 
with 7 = 500, L — 650 km and 50 kton of totally active scintillator detector. 
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Detection 
Efficiency 


Years 




Neutrino Beam (NH true) 


Anti-neutrino Beam (IH true) 


90% C.L. 


3a C.L. 


90% C.L. 


3a C.L. 


60% 


Ft A/P^T*^ 


250 
350 

(J (J V 

500 


8.1 x 10~ 3 

1.6 x 10~ 3 

6.7 x 10" 4 


2.8 x 10~ 2 
5.2 x 10~ 3 
2.1 x 10~ 3 


8.3 x 10~ 3 
1.3 x 10~ 3 
4.5 x 10" 4 


3.1 x 10~ 2 
4.3 x 10~ 3 
1.5 x 10~ 3 


1 f) vpars 


250 
350 
500 


4.2 x 10~ 3 
8.7 x 10~ 4 
3.7 x 10" 4 


1.4 x 10~ 2 
2.7 x 10~ 3 
1.1 x 10~ 3 


4.1 x 10~ 3 
6.7 x 10" 4 
2.4 x 10~ 4 


1.4 x 10~ 2 
2.2 x 10- 3 
7.4 x 10~ 4 


80% 


5 years 


250 
350 
500 


6.3 x 10~ 3 
1.3 x 10~ 3 
5.3 x 10~ 4 


2.1 x lO" 2 
4.0 x 10~ 3 
1.6 x 10~ 3 


6.2 x 10~ 3 
9.8 x 10~ 4 
3.5 x 10~ 4 


2.2 x 10" 2 
3.2 x 10~ 3 
1.1 x 10- 3 


10 years 


250 
350 
500 


3.4 x 10~ 3 
6.9 x 10~ 4 
2.9 x 10~ 4 


1.0 x 10~ 2 

2.1 x 10~ 3 
8.4 x 10~ 4 


3.1 x 10~ 3 
5.1 x 10~ 4 
1.8 x 10~ 4 


1.1 x 10~ 2 

1.6 x 10~ 3 

5.7 x 10- 4 



Table 4: Upper limit on sin 2 2#i3 which can be imposed with 5 years and 10 years data respectively, 
at the 90% and 3a C.L., for different benchmark values of 7. Columns 3 and 4 show the sensitivity 
using the v e beam and assuming normal hierarchy as true and columns 5 and 6 are the sensitivity 
possible with the v e beam and assuming inverted hierarchy as true. Results are shown for 60% 
and 80% ^ detection efficiency. 



significantly better than what one can achieve with the CERN-INO /3-beam set-up [19, 21]. 

We reiterate that the tremendous sin 2 29 13 sensitivity of the CERN-INO /5-beam set-up comes 
not really from high /5-beam flux, but from the fact that one has near-resonant matter effects 
at the near-magic baseline. While the latter feature smothers the problems arising from clone 
solutions, the former lends the experiment an extra sensitivity due to the dependence of matter 
effects on sin 2 2^13. This is why the sin 2 2^13 sensitivity depends on the true hierarchy. We have 
checked that if we run the experiment for 10 years in the neutrino mode alone and with inverted 
hierarchy assumed as true, then the sensitivity we expect is a meagre sin 2 2^13 > 1.8 x 10~ 2 at 
the 90% C.L. Therefore, a hierarchy independent sin 2 2#i 3 sensitivity can in principle be obtained 
after runs in both the neutrino and the antineutrino mode, unless one is very lucky. We will return 
to this issue again in section 8. 

6.2 Precision in measurement of sin 2#i3 

In this subsection we assume that sin 2 2# 13 (true) is large enough such that, for a given boost 7, we 
have a positive signal in the ICAL detector. For each value of sin 2 2# 13 (true), varied between the 
limiting value obtained in the previous subsection and the current upper limit, we generate the 
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Figure 9: Precision in the measurement of sin 2 2#i3 expected in the CERN-INO /3-beam experiment 
with 5 years running in the neutrino channel assuming normal hierarchy. Upper panels show the 
90% and 3a band of "measured" values of sin 2 2^13 for every sin 2 2^13 (true). Lower panels show 
the corresponding 90% and 3a CL. precision as defined in the text, as a function of sin 2 2^13 (true). 
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Figure 10: Same as Fig. 9 but for 10 years of running. 
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data and fit it back using the x 2 analysis to give us the range of "measured" values of sin 2 20i 3 . 
We plot in the upper panels of Figs. 9 and 10, the 90% and 3a C.L. bands of sin 2 20 i3 as a function 
of sin 2 2^13 (true) for 5 years and 10 years run respectively. The left, middle and right-hand panels 
are for 7 = 250, 350 and 500, respectively. Plots are shown for the v e /5-beam assuming normal 
hierarchy as true. In order to quantify the precision which can be obtained on sin 2 20i3, we define, 
for any value of sin 2 2^13 (true), 

„ sin 2 20i 3 (max) — sin 2 20iq (min) „ 

Precision = , ' 3 . 2 ' 3 , x 100 % , 27 

sin 20i3 (max) + sm 2(Ji 3 (mtnj 

where sin 2 29i 3 (max) and sin 2 26 13 (min) are respectively the largest and smallest values of sin 2 20 i3 
allowed at the given C.L. We show in the lower panels of Figs. 9 and 10, the variation of this 
precision function with sin 2 20 13 (true) for the 90% and 3a C.L. for 5 years and 10 years of running 
respectively, in the neutrino channel assuming normal hierarchy. Similar results are expected if 
we run the experiment in the antineutrino mode for inverted hierarchy. 

As expected, the precision in the sin 2 20i 3 measurement improves as sin 2 20i 3 (true) increases. 
However, a very interesting feature in the plots of Fig. 9 is the emergence of minima in the pre- 
cision, thereby implying that for values of sin 2 29 13 (true) beyond this limiting value, the precision 
actually deteriorates. 



7 Sensitivity to Sgn{/\m\^) 

We next turn our attention to the determination of the neutrino mass hierarchy. If the true value 
of 0i3 has indeed been chosen to be large by Nature, we can expect sizable matter effects in this 
experimental set-up, giving us a handle on the mass hierarchy. In our analysis, we assume a certain 
true mass hierarchy and generate the prospective data in INO-ICAL for the benchmark oscillation 
parameter values given in Table 2 and for different values of sin 2 20i 3 (true). We next fit this data 
with the wrong hierarchy and calculate a x 2 1 after marginalizing over (Am^l, sin 2 2 2 3 , 5cp as 
well as sin 2 20i3, as described in section 5. The results of this analysis are presented in Fig. 11. 
The figure shows, as a function of 7, the minimum value of sin 2 20 i3 (true) for which the wrong 
hierarchy can be disfavored at the 90% C.L. (solid curve for 5 years running and dotted curve for 
10 years running) and 3a C.L. (long dashed curve for 5 years running and dashed-dotted curve for 
10 years running). The left-hand panel is for a v e /5-beam and assuming that normal hierarchy is 
the true hierarchy, while the right-hand panel is for an v e /5-beam with the assumption that the 
inverted hierarchy is true. Of course, since large matter effects appear in only either the neutrino 
or the antineutrino channel for a given hierarchy, the sensitivity of the experiment to sgn(Am^ 1 ) 
would be zero if we used a u e beam when the true hierarchy was normal or an u e beam when the 
true hierarchy was inverted. 

In Table 5 we show the limiting values of sin 2 20i3(true), for which the wrong hierarchy can be 
disfavored at the 90% and 3a C.L. for three specific values of 7. We show the results for 5 years and 
10 years of running of the experiment in the neutrino and the antineutrino mode. Table 5 shows 
the results for both 60% and 80% ^ detection efficiency This can be directly compared with 
the expected potential of the other planned/proposed long baseline experiments in discriminating 
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Figure 11: The minimum value of sin 2 2^13 (true) as a function of the boost factor 7 at which the 
wrong hierarchy can be disfavored at the 90% and 3a CL. Results are presented for 5 years and 10 
years of running in the u e (left-hand panel) and v e (right-hand panel) mode, if the true hierarchy 
was normal and inverted respectively. 

between the right and the wrong hierarchy. The T2K experiment [53] has a baseline of only 295 km 
which entails almost negligible matter effects. The NOz/A experiment [54] on its own might find 
it hard to say anything conclusive about the neutrino mass hierarchy, unless Nature conspires to 
arrange oscillation parameters in a conducive way. The so-called T2KK experimental set-up with 
the 4 MW beam from the JPARC facility and two detectors, one in Japan and another in Korea 
[35] 20 , is expected to rule out the wrong hierarchy at the 3a CL. if sin 2 2#i 3 (true) > 5.5 x 10~ 2 for 
a total of 8 years of running, 4 years in the v e and 4 years in the v e mode. The CERN-MEMPHYS 
/3-beam and SPL superbeam proposals [42] again concern a baseline of only 130 km, which is too 
short to enable these beam experiments to decipher the mass hierarchy on their own. A higher 7 
/3-beam from CERN to (say) Gran Sasso, with a baseline of 732 km could have sensitivity to the 
mass hierarchy [24, 57]. It was shown in [56] that for 7 = 200 and L = 520 km, the wrong hierarchy 
could be disfavored at 3a CL. if sin 2 2^13 (true) > 2 x 10~ 2 after 8 years of simultaneous running 
in the neutrino and antineutrino mode and with a 500 kton water detector. Therefore we see that 
our sensitivity to the neutrino mass hierarchy can be well surpassed only by a Neutrino Factory, 
where the right normal hierarchy could be determined at 3a CL. for sin 2 2# 13 (true) > 3 x 10~ 4 , if 
the experiment is done at the magic baseline L = 7500 km. 

In Fig. 12 we show the minimum sin 2 2#i 3 (true) as a function of the running time T (in years) 
needed in this experimental set-up to discriminate between the wrong and the right hierarchy at the 
90% and 3a CL. The left-hand panel is for a run in the v e mode with true normal hierarchy, while 

20 For an alternate T2KK proposal, see [36]. 
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Detection 
Efficiency 


Years 




Neutrino Beam (NH true) 


Anti-neutrino Beam (IH true) 


90% C.L. 


3a C.L. 


90% C.L. 


3a C.L. 


60% 


Ft A/P^T*^ 


250 
350 

(J (J V 

500 


3.7 x 10~ 2 
9.7 x 10~ 3 
7.7 x 10~ 3 


6.6 x 10~ 2 
1.5 x 10~ 2 
1.0 x 10~ 2 


4.7 x 10" 2 
9.2 x 10~ 3 
5.9 x 10~ 3 


8.7 x 10~ 2 
1.4 x 10" 2 
8.2 x 10~ 3 


1 f) vpars 


250 
350 
500 


3.0 x 10~ 2 
8.3 x 10- 3 
6.6 x 10~ 3 


4.7 x 10~ 2 
1.2 x 10~ 2 
8.7 x 10~ 3 


3.8 x 10~ 2 
8.1 x 10~ 3 
5.4 x 10~ 3 


6.0 x 10~ 2 

1.1 x 10- 2 
7.0 x 10~ 3 


80% 


5 years 


250 
350 
500 


3.4 x lO" 2 

9.3 x 10~ 3 

7.4 x 10- 3 


5.7 x 10- 2 
1.3 x 10~ 2 

9.8 x 10~ 3 


4.3 x 10" 2 
8.9 x 10~ 3 
5.9 x 10~ 3 


7.3 x 10- 2 
1.3 x 10~ 2 
8.2 x 10- 3 


10 years 


250 
350 
500 


2.8 x 10~ 2 
8.1 x 10~ 3 
6.6 x 10~ 3 


4.2 x 10~ 2 
1.1 x 10~ 2 
8.5 x 10~ 3 


3.5 x 10~ 2 
7.8 x 10~ 3 
5.4 x 10~ 3 


5.4 x 10" 2 
1.0 x 10~ 2 
6.9 x 10" 3 



Table 5: Minimum value of sin 2 2^13 (true) for which a CERN-INO beta-beam experiment can rule 
out the wrong hierarchy at the 90% and 3a C.L. after 5 years and 10 years of running, for different 
benchmark values of 7. Columns 3 and 4 show the sensitivity using the v e beam assuming normal 
hierarchy as true and columns 5 and 6 are for the v e beam assuming the inverted hierarchy as 
true. Results are shown for 60% and 80% ^ detection efficiency. 



the right-hand panel shows the results for running in the v e mode with true inverted hierarchy. 
We show the results for the three chosen values of 7. We see, for example, that if sin 2 2# 13 (true) = 
5 x 10~ 2 and the normal hierarchy was true, then the wrong inverted hierarchy could be disfavored 
at the 3a C.L. after just 9 (26) months of running in the v e mode with 7 = 500 (350). 

Fig. 13 shows the impact of 8c p on the hierarchy sensitivity of this experiment. The left-hand 
panel is for the neutrino channel with normal hierarchy as true, while the right-hand panel is for 
the antineutrino channel with inverted hierarchy as true. We show the 90% (dashed curves) and 
3a (solid curves) C.L. curves for 5 years (thick lines) and 10 years (thin lines) running of the 
/3-beam with 7 = 350. We note that loss in hierarchy sensitivity due to the uncertainty in 5cp 
is very marginal, owing to the fact that we are proposing this set-up at a near-magic baseline. 
The best sensitivity to hierarchy comes for 5cp — 125 in the neutrino mode and 5cp — 300 in the 
antineutrino mode. 

8 Running the v e and v e /3-Beam Simultaneously 

An exquisite feature of magnetized iron detectors is their charge identification capability. Thus, 
if both the v e and the v e /5-beam were produced simultaneously at CERN, the oscillated could 
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Figure 12: The minimum value of sin 2 2^13 (true) as a function of the running time necessary to 
disfavor the wrong hierarchy at the 90% and 3cr C.L. in the neutrino mode assuming normal 
hierarchy as true (left-hand panel) and in the antineutrino mode assuming inverted hierarchy as 
true (right-hand panel). Results are shown for three values of 7. 



be measured distinctly from the oscillated at INO 21 . Since near-resonant matter effects appear 
only in either the neutrino or the antineutrino channel, this gives the experiment an added edge 
to unravel both the mass hierarchy as well as #13 simultaneously, within a much shorter time. The 
argument goes as follows. 

Throughout the paper so far, we had been tacitly working in a scenario where we have just 
one type of /3-beam from CERN and the hierarchy was conducive to having large matter effect, 
the key behind the excellent #13 as well as mass hierarchy sensitivity. However, we have to bear in 
mind that in the absence of any prior knowledge of the true value of #13 and/or the mass hierarchy, 
one would not know which of the two modes - neutrino or the antineutrino - would be preferable 
for extracting maximal matter effects. With luck, it might so happen that the beam is run in 
the neutrino mode first and the true hierarchy turns out to be normal and sin 2 2# 13 (true) not 
too small. In that case, we would see large number of events in the detector and know without 
ambiguity both the mass hierarchy as well as the value of #13. However, if the experiment is run 
in say only, the neutrino mode and if we do not see enough number of events (as expected from 
large matter effects), it would be difficult to tell if the non-observation should be attributed to 
the true hierarchy being inverted or whether the mixing angle #13 is very small. We would find 
ourselves in a similar impasse if a v e /3-beam is used and one observed a very small number of 

21 The v e and i> e /3-beam fluxes could be produced in distinct bunches from the storage ring. Therefore, in 
principle, it should be straight-forward to distinguish the neutrino from the antineutrino events at any detector 
using the nano-second time resolution that most detectors possess. However, since ICAL will be magnetized, it 
could directly distinguish between the two channels irrespective of the timing information. 
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Figure 13: The effect of Sep uncertainty on the limiting value of sin 2 2^13 (true) for which the 
hierarchy determination will be feasible for 7 = 350. For all values to the right of the solid curves 
(dashed curves) it would be possible to rule out the wrong hierarchy at the 3a (90%) C.L. Thick 
(thin) lines are for 5 (10) years of running. Left-hand panel is for v e /3-beam assuming normal 
hierarchy is true, while the right-hand panel is for v e /3-beam assuming true inverted hierarchy. 

events. This deadlock would then be broken only if we run the /3-beam in both the neutrino and 
the antineutrino mode. 

Recall that a /3-beam facility has the useful feature that it can be operated simultaneously in 
the v e as well as v e mode, since both the S B and 8 Li ions are positively charged, and therefore 
move in the same direction in the storage ring. The only constraint is that, since the magnetic 
field would be the same for both channels, the different charge-to-mass ratio for the two candidate 
ions would fix the boost for the antineutrino channel to be 0.6 times that in the neutrino channel. 
Therefore, we could work within a set-up where we have simultaneous v e and v e beams with say 
7 = 500 and 300 respectively for the two channels. Since INO-ICAL would have excellent charge 
identification capability, it would distinguish the /i~ produced from oscillated and ji + produced 
from oscillated u^, when the v e and v e beams are run simultaneously In this scenario, unless 
#13 is extremely small, we would observe an enhancement in the number of events due to earth 
matter effects in one of the channels, irrespective of the true mass hierarchy. Therefore, a direct 
comparison of the number of fi~ and fi + events would immediately tell us the mass hierarchy for 
#13 not small. If the true value of #13 is indeed chosen to be very small by Nature, then we would 
see no enhancement in the event rate in either the neutrino or the antineutrino channel and one 
would be able to place an upper limit on the value of sin 2 2#i 3 unambiguously. 

We have checked that with v e and v e beams running simultaneously with 7 = 500 and 300 
respectively for the two channels, we could obtain a #13 sensitivity of sin 2 2#i 3 < 4.5 x 10 -3 at 3cr, 
irrespective of the true mass hierarchy, after 10 years of running. We would have sensitivity to 
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the neutrino mass hierarchy for at least sin 2 2#i3 > 2.1 x 10~ 2 at 3a after 10 years of simultaneous 
running, again irrespective of whether the true hierarchy was normal or inverted. One could argue 
that running the beam successively with 7 = 500, first in the v e mode followed by the V e mode 
would give an unambiguous #13 sensitivity limit of sin 2 20 13 < 1.5 x 10~ 3 at 3a and the mass 
hierarchy sensitivity limit of sin 2 26i% > 8.2 x 10~ 3 at 3a. However, note that it is not possible 
to run the v e beam with the same 7 as the v e beam without any accelerator upgrade. Running 
the /5-beam simultaneously in both polarities has the additional advantage that if sin 2 2# 13 (true) 
is indeed large in Nature and within reach of the CERN-INO /9-beam set-up, then one would have 
a surplus of either fx~ or /i + events in INO-ICAL. This would directly tell us the neutrino mass 
hierarchy unambiguously. 

9 Conclusions 

We explored the physics potential of an experimental program with the /3-beam source [y e or u e ) 
at CERN and a 50 kton magnetized iron detector at INO. The CERN-INO distance is very close 
to the magic baseline, which is ideal for performing a degeneracy-free measurement of the mixing 
angle 13 and the neutrino mass hierarchy. The large baseline also opens up the possibility for large 
matter effects. For the CERN-INO baseline, the largest earth effects corresponding to resonant 
conversion of (anti) neutrinos in matter, happen for (anti)neutrino energy of E ~ 6 GeV. or 

produced at these energies through oscillations can be effectively detected in magnetized iron 
calorimeters via the ^ they produce, the detection energy threshold being about 1 GeV. For the 
most common ions used in /5-beam studies, this range of (anti) neutrino energies can be achieved 
for Lorentz factor 7 ^ 10 3 . In this paper we have considered the new possibility of using the S B 
and 8 Li as source ions for v e and v e /3-beam respectively. Owing to the large end-point energy of 
these ions, the resultant energy spectrum of the Lorentz boosted /9-beam for these ions is harder 
by more than a factor of 3, compared to the standard alternative ions 18 Ne and 6 He. Therefore, 
with these ions we would get a /5-beam flux peaked at about 6 GeV, for plausible values of the 
Lorentz boost 7 in the range 250-500. We showed that for this /3-beam flux from CERN, it will 
be possible to get an essentially background free measurement of the "golden channel" probability 
P eM , through observation of fi~ (or /i + ) events in the ICAL detector at INO. We argued that the 
number of events in this experimental set-up is dictated directly by the extent of near-resonant 
earth matter effect, which significantly enhances the transition probability P efl . The extent of 
earth matter effect in turn is governed by the neutrino mass hierarchy and the value of sin 2 2#i 3 . 
The CERN-INO /3-beam experiment that we propose therefore emerges as a powerful tool to pin 
down the neutrino mass hierarchy and 9 13 . 

We simulated the prospective data in this experimental set-up for 5 years and 10 years of 
running of the /3-beam in either the neutrino or the antineutrino mode and presented the sensitivity 
results from a rigorous x 2 analysis, after marginalizing over all the oscillation parameters and 
taking into account the systematic uncertainties coming from the /3-beam source, the detector, 
as well as the theoretical calculation of the interaction cross sections. For 7=500, if true value 
of $i3 was zero, this experiment with 60% detection efficiency and 5 years of running would put 
a limit of sin 2 26 13 < 2.1 x 10~ 3 (sin 2 26 13 < 1.5 x 10~ 3 ) at 3a C.L. with the v e (p e ) /3-beam for 
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the normal (inverted) mass hierarchy. With 80% detection efficiency and 10 years of running the 
corresponding limits would be sin 2 29 u < 8.4 x 10" 4 (sin 2 26 13 < 5.7 x 10~ 4 ) at 3a C.L. with the v e 
(p e ) /5-beam for the normal (inverted) mass hierarchy. With 60% detection efficiency and 5 years 
of running the wrong hierarchy could be ruled out at the 3a C.L. if sin 2 2# 13 (true) > 1.0 x 10~ 2 
(sin 2 2#i 3 (true) > 8.2 x 10 -3 ) with the v e (p e ) /5-beam for the true normal (inverted) mass hierarchy. 
If the detection efficiency was 80% and running period 10 years the corresponding sensitivity would 
be enhanced to sin 2# 13 (true) > 8.5 x 10~ 3 (sin 2 2# 13 (true) > 6.9 x 10~ 3 ) respectively at the 3a 
C.L. We showed that if sin 2 20 13 (true) = 0.05, 7 = 500 and the true hierarchy was normal, then the 
inverted hierarchy could be disfavored at 3a C.L. with a v e /9-beam within just 9 months of running. 
Finally, we considered the scenario where the v e and v e flux could be produced simultaneously at 
the CERN /3-beam facility. Since ICAL at INO would be magnetized and hence would possess 
charge identification capability, we argued that an unambiguous information on #13 and the mass 
hierarchy could be obtained within a much shorter time scale if sin 2 26*i 3 (true) is large. If #13 is 
very small, one would not be able to determine the mass hierarchy, but could still save time in 
putting an upper bound on the mixing angle #13. 

In conclusion, the sensitivity for sin 2 26>i 3 and the neutrino mass hierarchy that is expected 
in the CERN-INO /9-beam experiment appears to be better than most of the other experimental 
proposals involving superbeams and /9-beam and can only be surpassed by the NuFact-II experi- 
ment. 
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